On the Origin of Quantum Oscillations in the Mixed State of Anisotropic 

Superconductor 
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For the vortex lattice in an anisotropic superconductor 
with well-separated cores (H c i <C B <C H C 2) it is shown that 
sizeable de Haas- van Alphen oscillations are caused by the lev- 
els' crossing of the energy threshold separating localized and 
extended states of excitations moving in the average magnetic 
field, B. 

PACS numbers: 74.20. Fg, 72.15. Gd, 74.60. -w 

Recent experiments [1,2] show that the de Haas- van 
Alphen (dHvA) effect persists in the superconducting 
(SC) state for the magnetic field, B, as low as B ~ 
0.3 ~ 0AH c2 . The effective SC-gap, A, at such fields 
is large enough to preclude motion of an electron along a 
closed Larmour orbit with some radius, rr, ~ vp/ui c , ex- 
ceedingly larger than the coherence length, £ , and the in- 
tervortex distance, d [3]. The dHvA-signal is expected to 
weaken exponentially as exp(— A/uj c ) [4]. In that which 
follows we suggest a new mechanism for the quantum 
oscillations in the mixed SC-state. 

In normal state the dHvA-effect is brought about by 
levels crossing the chemical potential, fx, with the field 
variation. The oscillations are periodic in B~ x because a 
minor field change, AB/B ~ u) c J /i, is enough to push a 
level across \i. 

Electron- or hole-like character of SC-excitations de- 
pends on the extent their energy exceeds the gap. Even 
for a "d-wave" superconductor [5] levels cannot cross the 
chemical potential. In this sense, there is no difference 
between a "d-wave" or any other anisotropic SC. 

It is shown below that a new energy threshold takes 
over the role of the chemical potential in the SC-state. 
Consider, for example, an anisotropic superconductor 
with a spectrum, e(p) = y 1 v 2 F (j> — pf) 2 + | A(p) | 2 . As- 
sume, for simplicity, that the gap, A(p), has only one 
maximum, A max , and one minimum, A min , along the 
Fermi Surface (FS). Excitations with e(p) > A max have 
itinerant behavior, while at A max > e(p) > A m i n this 
is only true for excitations with a proper p. The lat- 
ter become localized in a magnetic field, for the Lorentz 
force changes the p-direction. Excitations with the en- 
ergy larger than A max , may move along an extended Lar- 
mour orbit. 

To pose this phenomenon as a theoretical prob- 
lem, consider the limit of well separated vortices d 
£o(Hd -C B <C H c2 ). The vortex cores occupying only 
a minor fraction of the volume may be neglected. A 



typical electron trajectory would run across the "bulk" , 
where the gap amplitude is saturated: 



A(p,r) ^ A(p)exp(^(r)) 



(1) 



The method [5] to treat the problem is based on averaging 
the Gor'kov system of Green functions over quasiclassical 
trajectories [6] (all notations below from ref. [5]). The 
Green functions are presented in terms of the position, 
tp, of an electron along the FS. The core of the method 
is given by eqs. (23-28) [5]. The Gor'kov matrix being 
diagonalized, the whole problem reduces to solving the 
following Schrodinger equation: 



+ \A 2 ^) - u c A'(tp)]y = E 2 y 



(2) 



The term, — u> c A'(tp), plays no role and will be omitted. 
In (2) we have also left out the term h(tp) of eq. (31) [5]. 
The Doppler shift (31) [5] is essential for the magnitude 
of the effect, and will be taken into account later. 

The eigen functions are given by solutions of (2) satis- 
fying the periodicity condition (tp — > tp + 2tt) for 



y(tp)e-^ 



(3) 



where n = n/u) c comes from presenting the chemical po- 
tential in the form [5] : 



/z = lu c N + fi 



(4) 



With fi being large, Nq 3> 1, and the specific Nq does not 
affect the pattern of periodic (in B~ r ) oscillations of the 
magnetization which depends on k in the interval (0,1). 
Eq. (2) and (3) become the problem of finding the band 
structure for a particle moving in the periodic potential 
A 2 (tp) with k as a quasimomentum. 
Re-write (2) in the form: 



(A 2 M 



)y = (e 2 



(2') 



and consider (2') first in the quasiclassical (WKB)- 
approximation (oj c « A). At \E\ ^> A the periodicity of 
(3) leads to the spectrum of free electrons in the magnetic 
field: E n — uj c n + fj. At E 2 — A 2 nax < the attractive 
potential in (2') has many (A/u> c 3> 1) "localized" levels, 
(tunneling across the barrier is neglected). The bound- 
ary separating "extended" and "localized" states in the 
WKB-sense lies at A ma . x . Introduce in (2'): 



E 1 



(5) 



l 



for \E\ close to A max . The WKB-solutions are [7]: 

y ± (<p) = AiS'^)) 1 ' 2 cxp[±z f S'(<p)d<p] (6) 
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with 

UcS'ifp) = [2A max (-e) + Al ax A 2 (p)} 1/2 (7) 

(A is the normalization coefficient). The BCS-factors, 
u(ip) and v((f), in (23) [5] are to be normalized together: 
\u 2 \ + \v 2 \ = 1, (the bar in (. . .) means the normalization 
integral: (2-k)^ 1 J" * (. . .)d(^). Two auxiliary expressions 
which follow from eqs. (26-28) [5]: 

W\ = (1/2) {W+ (tuJ c /2E)(y*y> -yy*')} 

W\ = (1/2) { W - ^ c /2i?)(r?/-yr')} (8) 

immediately show that \y\ 2 = 1. 

Expression (6) [5], containing oscillatory effects 

M = -^E» ( 9 ) 

A 

(A enumerates the eigenvalues, factor 2 added to (6) [5], 
accounts for spins, and n(E\) = 1 at E\ < and T = 0), 
becomes an integral over A with the use of the Poisson 
formula: 

+oc +oo 

£f(A-»)= 2 (10) 

— CO /c— — CO 

Integration by parts transforms M osc into [5]: 

- /T^s (^ dA <101 

Although integration over A acquires the meaning only 
after the connection between A and the energy is es- 
tablished, the threshold separating "localized" and "ex- 
tended" (in the WKB-sense (6,7)) states is already seen 
in eq. (10'): for "localized" states, \E\ < A max , the 
wave functions are real, and from (8) |u 2 | = 1/2. For 
"extended" states (6) we have: 

KF=i|i±|i|A A | 2 } (ii) 

with I^Ap-a function of energy (see (6)). The derivative 
in (10') thus eliminates states below Amax 

with \u\\ 2 = 

1/2 being energy independent. 

Returning to summation over A, eqs. (10) and (10'), 
we need to construct such the function: 

X(E) = $(£)/2tt (12) 

that the provision: 



$(£„) = 2nn (12') 

would enumerate all energy levels in the consecutive or- 
der. In the WKB-approach &(E) is given by: 

S(2w, -e) = (l/w c ) / S'fa -e)dp - 2-kk (13) 
Jo 

which at large energies matches the Landau free electron 
spectrum. The approach falls short near (— e) = 0. 
Choose A 2 (<p) near A max as 

AV) = AL*(1-<V) (14) 
Expanding (13) in (— e) > 0, one obtains: 

5(271-, -e) ~ 5(2tt, 0) - (1/2) In(//A) (15) 
with the useful notation in (14): 

I = {-2e)/uj c a 1/2 ■ A = (a 1 ' 2 A max /oj c ) > 1 (16) 

Similarly, the factor \A\ 2 in (6) and (11) at small I is 
proportional to: 

\A\ 2 cx [HI/A)}- 1 (15') 

Because of the singularity (15) the "numbering" func- 
tion $(£?) cannot be comprised of the two WKB- 
branches, the one that is given by (13) (at (— s) > 0), and 
the other which counts "localized" states ((— e) > 0). 

Note that far away from ip = the WKB-solution 

y(ip) = ay+ (<p) + by^ (ip) (17) 

is still correct. With the use of (5,14,16), eq. (2') can 
be solved near ip = in terms of the parabolic cylinder 
functions. It establishes the matrix relation between co- 
efficients (a, b) in (17) on the R.H.S. of <p — + , and the 
other set (a',b') on its L.H.S., ip = 0_: 




Beginning at ip = 0+, moving along with eq. (17) toward 
(27r)_ and using (18), the periodicity condition for (3) 
provides the equation: 

R(l) = \a\(e i§ + e~ i§ )= 2 cos2ttk (19) 

Abbreviations in eq. (19) are: 

S = S(2tt, -s) -6, a= \a\ exp(i0) (19') 

Of the two solutions of eq. (19) we chose 

e lS = j^j {cos 2tth + i^\a\ 2 - cos 2 2ttk} = p(l) (20) 
\a\ 

because at (— e), i.e. I, large and positive \a\ => 1, 0, 
and we return in this limit to (13). On the real axis of I 
the function 



= S(2n,l) 



0(1) - -lnp(i) 



(21) 



is positive, with d&/dl > 0, and matches asymptotically 
(at I — ► oo ) the free electron spectrum. With the help of 
(12) and (21) the oscillatory magnetization is 



M osr . = 



~2^ 2 ~c 



dl (22) 



$(/) being continued analytically into the complex l- 
plane, the integration may be shifted into the upper half- 
plane (from (11, 15') one has the behavior of d(\ui\ 2 )/dl 
at large \l\). Consider singularities in (22). Thus from 
the expression for a| [8]: 



(1 



-ttZn 1/2 



we conclude that branching points in (23) lie at 
l m = ±(2m+ l)i 



(23) 



(24) 



This is also true for 0(1) (See S(2n, I) below). The defini- 
tion of p(l) together with (23) for |a|, leads to the square 
root singularities at 



C = ±(2m- 



l)i ln(sin 2 2nn) 

TT 



Ira + l (24') 



For $(/) to be analytical in a strip at the real axis, the 
branch-cuts in the /-plane, caused by singularities (24, 
24'), must be chosen parallel to the imaginary axis. 

The integral in (22) may be bent to the contours, C\ 
and C2, each encircling the branch-cuts (24) and (24') 
in the upper half-plane. The non-analytical terms of eq. 
(15) at |A| ~ 1 are now absent in S(2n, l)(c ~ 1): 



S(2nJ) ~ S , (2tt,0) + ^[ln(A • c)] 

iL 



In 

2i 



r(; 



^2 - f ) 



(25) 



Both integrals (along C\ and C2) rapidly converge. 

It is necessary to normalize \ui(ip)\ 2 with the accuracy 
better than that given by the WKB-approximation in (6), 
(11). Fortunately the properties of the Bloch functions 
in a one-dimensional periodic potential are well-studied. 
With the help of eq. (4.18) in ref. [9] and our eqs. (8) 
we derive: 



|w;(<^)| 2 = - — na 1 / 2 sin 2ttk 



dR 

~dT 



(26) 



(Here R(l) is the R.H.S. of (19)). After differentiation on 
I in (22) had eliminated (1/2) in (26), one may use for 
the rest the rapid convergence of the integrals to integrate 
back by parts in (22). Single terms under the sum symbol 
(22) become: 



i-ira 1 / 2 sin27TK 



LK 



s 



exp(iK^(l))dl 
(sin 2 2ttk + e -Ti)i/2 



(27) 



with integrals running along C\, Ci- Assume that In A 3> 
1 in (25). First, the term in (22) with K = 1 prevails. In 
addition, as seen from (25), it is enough to consider the 
nearest singularities with m = in (24, 24'). Defining 
the branches of the square roots in (27) properly, from 
(27) one obtains contributions to I\ from the two contour 
integrals, over C\ and C2, correspondingly: 



h(Cx 



(2^f/ 2 a l 'Hg2 



(AcOVSQnAc) 1 / 



7rK e *S(2 w ,0)+i2 



(28') 



/i(C 2 ) 



ita 1 / 2 sin 2-KKe" 
(AcjV^hAc) 1 / 2 



r(i - %) 
r(i + %) 



1/2 



-e iS(2Tr,0) + ^ ln(Ac) 



(28" 



with l from (24'). (h has the form (28', 28") at l 
not too small (Iq ~ (InAc) -1 ). Otherwise the two con- 
tours C\ and C2 start to merge. Also, if l becomes large 
(sin 2 2nn — ► 0), expression (25) for S(2n, I) being correct 
at |Z| ~ 1, ceases to be applicable). 

Expressions (28) explicitly present the periodic (in 
B^ 1 ) oscillations in the magnetization as k varies in the 
interval (0,1). It is notable that the amplitude is of the 
order of (wc/Amm) 1 / 2 , i.e., is not exponentially small. 
Both (28') and (28") lead to the large content of the 
higher harmonics. In principle, M osc could be measured 
directly as a function of small changes in B. Ii(C\) dis- 
closes a rather regular behavior in AB (i.e., k), while 
-MC2) rapidly becomes chaotic due to the phase factor in 
(28"), exp [i(lo/2) ln(Ac)], contributing into higher har- 
monics. At the Fourier analysis of the dHvA-signal a few 
first harmonics are expected to be seen with the intensity 
of the order of 



(uJA^) 1 ' 2 



(29) 



Unfortunately, (29) does not take into account scatter- 
ing of electrons on the flux lines. The term h(tp) of eq. 
(31) [5], if included, adds to the potential of eq. (2'): 



2A r 



c h(<p) 



(30) 



Even though h(ip) ~ vp/d is small compared with A max , 
(30) drastically distorts the potential near ip = 0. It 
is a local h max in the vicinity of the maximum in A(ip) 
which now sets in the energy threshold between "local- 
ized" and "extended" states. Note that although h((p) is 
rather irregular (for a given trajectory) and does change 
typically on the scale of Sip ~ (duj c /vp), its local maxima 
produce potential barriers in (30) which remain impen- 
etrable in the quasiclassical sense. The above analysis 
of M osc can be performed in exactly the same manner 
as above around h max . There is a change in the scale 
(29), because the curvature, a, near a maximum in (30) 
is much higher than in (14). Without going into details, 
we comment that this only increases the effect, because 
the potential h(tp) nG&r hmax comprises 3. much sharper 
barrier, if compared with eq. (14). 



The major destroying effect comes from the phase 
factor in (25), S(2n,0). At h{<p) ^ 0, S(2tt, 0) may 
be expanded in 8h{ip) = h(ip) — h max . (Now (— e) => 
E - A max - hmax)- The fluctuating part, 6S(2n,0), is 

6S(2n,0) = -(A max /uj c ) f 6h(<p) {A 2 max - A 2 ^)Y 1,2 d v 

Jo 

Since (Sh((p)) (the average over all trajectories) is obvi- 
ously zero, fluctuations in exp(i5S(2Tr, 0)) lead, as in [5], 
to an effective Dingle factor of the form: 

exp(-v F /duj c ) exp(-(A/w c )(f /d)) ( 32 ) 

The exponent in (32) provides much more favorable con- 
ditions for observation of the dHvA-effect than previous 
results [4]. 

To conclude, in the developed mixed state of an 
anisotropic superconductor there exists the energy 
threshold sorting excitations into two categories: the lo- 
calized and extended ones. Crossing this threshold by 
the excitations' levels at the change of the magnetic 
field comprises the new mechanism for quantum oscil- 
lations. Scattering on the flux lines reduces the dHvA- 
effect. Nevertheless, the effect remans bigger than antic- 
ipated. 
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